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The analytical solutions for the general-four-wave-mixing hyperpolarizabilities \ ( — (loi +L02 + 
LU3); u>i, U2, ^3) on infinite chains under both Su-Shrieffer-Heeger and Takayama-Lin-Liu-Maki mod- 
els of trans-polyacetylene are obtained through the scheme of dipole-dipole correlation. Ana- 
lytical expressions of DC Kerr effect x' 3 ' ( — 0, 0, lu), DC-induced second harmonic generation 
x' 3 ' (— 2uj; 0, lj, lj), optical Kerr effect x' 3 ' ( — t*>; w, — U), w) and DC-electric-field-induced optical rec- 
tification x' 3 ' (0; to, — to, 0) are derived. By including or excluding Vk terms in the calculations, 
' comparisons show that the intraband contributions dominate the hyperpolarizabilities if they are 

included. Vfc term or intraband transition leads to the break of the overall permutation symme- 
try in x' 3 " 1 even for the low frequency and non-resonant regions. Hence it breaks the Kleinman 
symmetry that is directly based on the overall permutation symmetry. Our calculations pro- 
vide a clear understanding of the Kleinman symmetry breaks that are widely observed in many 
experiments. We also suggest a feasible experiment on x' 3 ' to test the validity of overall per- 
mutation symmetry and our theoretical prediction. Finally, our calculations show the follow- 
ing trends for the various third-order nonlinear optical processes in the low frequency and non- 
resonant region: X^_ rM (-3w; w . w . w ) > xion-r e3 (-2^; 0, w, to) > xil- res (-^^>-^^) > 
XnL-Tes(-u;0,0,u) > X™„_ re8 (0;w,-w,0), and in the resonant region: x^si-uj; 0, 0, to) > 
O [ XrV s (-to;to, -to,to) > Xr 3 l(~2aj;0,w, u) > X?es(Q; to; -to, 0) > 3w;w, u, to). 

PACS numbers: 78.66.Qn, 42. 65. An, 72. 20. Dp, 78.20.Bh 
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The nonlinear optical (NLO) properties of 7r-conjugated polymers have been extensively studied by both experimen- 
talists [1-11] and theorists [12-29]. Among these polymers, polyacetylene (PA) is the simplest conjugated polymer 
and has been extensively studied [30]. PA consists of chains of CH units that form a pseudo one-dimensional(lD) lat- 
tice. Classical periodic single electron models like Su-Shrieffer-Heeger (SSH) [31] and Takayama-Lin-Liu-Maki (TLM) 
[32] have been established to interpret the optical properties of polyacetylene. Recent experiments have measured 
the spectrum of third-order harmonic generations (THG) in polyacetylene [5,6], which requires theoretical explana- 
tions. The contribution from the lattice relaxation to the third-order optical susceptibilities (x ) is believed to be 
small because of the slow process involved in the lattice relaxation [33]. Based on either the periodic single electron 
models (SSH and TLM) or the inclusion of electron-electron interactions extended from both models, the frequency 
dependence of the THG in PA has been extensively studied theoretically [12-26]. 

Even for these simple single electron models, the theoretical studies on frequency dependence of x' 3 ' have been 
■ approached in various ways: (i). Based on the different gauges, this problem can be studied either from dipole-dipole 
(DD) [12-14,19-23] or from current-current (JJ) correlations [15,16]; (ii). Based on the size of chain, this problem 
can be treated either as a finite periodic chain [18-20,23-25] or an infinite one [12-15,21]. Furthermore, from the 
point of view of different techniques to utilize the perturbative method, this problem can be solved either by quantum 
mechanics [12,17-24] or by the field theory [13-16]. Obviously, different approaches on the same model should yield the 
same results if every step is carried out correctly. However, this seemingly simple problem turns out to be nontrivial 
in the actual computations because different results of x^ for polyacetylene have been obtained from the different 
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approaches. A typical example is the THG calculation. There different gauges yield different results both numerically 
and analytically [13-15,19]. On the other hand, in the formulation of DD correlations, there are some controversies 
about whether or not the Vk term should be included for an infinite periodic system [28,29,34-36]. It is our interest 
to illustrate the differences and understand the reasons why the discrepancies exist in the theoretical calculations. 

Recently, many experiments showed the general failure of Klcinman symmetry [37] in describing the low frequency 
and non-resonant nonlinear optical properties of many materials and practical systems [38]. But most experiments 
only measured the x^ m systems without centro- or inversion symmetry. Based on the possible mutual exclusion 
property between Krammcrs-Kronig dispersion relations and Klcinman symmetry, Dailcy, Burke and Simpson asserted 
the general failure of Kleinman symmetry. However, physically Kleinman symmetry is only a direct consequence of the 
overall permutation symmetry [37,39]. Thus we would also like to investigate the overall permutation symmetry under 
SSH and TLM models in this work. For the trans-polyacetylene system where the centro- or inversion symmetry is 
applied, x^ vanishes. Hence we will use the non-trivial results for x' 3 ^ to discuss the validity of both overall 
permutation and Kleinman symmetries. 

In this and the subsequent [40] papers, we will present the analytical results of hyperpolarizabilities (x' 3 ^) for both 
the SSH and TLM models based on the field theory. Both SSH and TLM models describe periodic single electron 
systems. The analytical form for the nonlinear optical response can be obtained under both models and used to 
illustrate the subtle issues in periodic systems. We will make detailed comparison between our calculation and other 
theoretical results of hyperpolarizabilities in both papers. The purpose of present paper is to elucidate the physical 
contribution of each term in the formulation of dipole-dipolc correlations. In these models, the important role played 
by the operator Vk due to its noncommuting feature with ordinary functions of k is evident. Our calculations show 
that it is the Vk term, which is physically related to the intraband transition, that leads to the break of the overall 
permutation symmetry and Klcinman symmetry in the non-resonant region. Based on our calculations, we also suggest 
some experiments on the measurement of x f° r certain physical systems to test the validity of overall permutation 
symmetry and our theoretical predictions. 

This paper is organized as follows. In Section II, necessary analytical tools and general theoretical framework are 
introduced. In Section III, analytical results of four- wave- mixing (FWM), DC Kerr effect (DCKerr), DC-induced 
second harmonic generation(DCSHG), optical Kerr effect or intensity-dependent index of refraction (IDIR), DC- 
electric-field-induced optical rectification (EFIOR) are derived. The discussions on those results and symmetries, 
comparisons with other theories and some suggested experiments are then followed in Section IV. A brief conclusion 
is presented in Section V. Finally, the details of analytical calculations are given in the Appendix. 



II. THEORY 



A. Models and dipole operator 

Both SSH (or dimerized Huckel models) and TLM models have been thoroughly studied in the irans-polyacetylene 
problems [30]. The SSH model is an infinite ID periodic model, and described by the following Hamiltonian [31]: 

Hssh = - E k + (- 1 ) t1 (All, A, + CjA+i,.), ( 2J ) 



l,s 



where to is the transfer integral between the nearest-neighbor sites, A is the gap parameter and Cj g (Cj, s ) cre- 
ates (annihilates) an 7r electron at site / with spin s. For the SSH model, each site is occupied by one electron. In the 
continuum limit, the SSH model will tend to the TLM model given by the formula [32]: 

Htlm = {ia 3 v F d x + Aoi) ¥(x), (2.2) 

where &(x) = (^{(x) , ^(x)) is the two-component spinor describing the left-going and right-going electrons, vf is 
the Fermi velocity, A is the gap parameter, and m (i — 1, 2, 3) are the Pauli matrices. 
Under the DD correlation, the interaction Hamiltonian is expressed by the formula 

H E . r = -eE f = — D E, 

with e the electron charge and E the electric field described as follows: 

E(r,t) =Eoe 4k - r -^ 4 , (2.3) 
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where Eo is the amplitude, k and u> are the wave vector and frequency, respectively. 

For periodic systems, the position operator r is often conveniently defined in the momentum space [44]: 



r„ k ,„'k' = ^<5„,„<V k (5(k - k') + Q n ,„,(k)<5(k - k'), (2.4) 

where 

0„ n /(k) = - / u* n k (r)V k u„' jk (r)dr, (2.5) 
v J v 

with v the unit cell volume, M njk (r) the periodic function under the translation of lattice vector [41]. Obviously, 
u Bj k(r) is related to the wavefunction ip of Bloch states by the formula: 

VVk(r) = u n , k (r)e^ r , (2.6) 

where n and k are the band index and crystal momentum, respectively. 

The energy bands of ID SSH model are simple - conduction and valence bands. Thus, this ID model avoids 
certain problems due to the discontinuous and non-analytical feature of Bloch wavefunctions for a composite band in 
higher-dimension periodic models [42] . 

Following the same procedures described in previous work [13,14,43], we consider the momentum space represen- 
tation of the Hamiltonian given by Eq.(2.1). With the aid of the spinor description s (t) = (a { k c s (t) , a]^ s (t)), where 

aj. c s (£) and (t) are excitations of electrons in the conduction band and the valence band with momentum k and 
spin s, we obtain the following formula: 

HssH(k,t) = H a + H E .r (2.7) 
= ]T e{k)^ ktS {t)a^ s {t)-D-E e^\ (2.8) 

where a are the Pauli matrices, the operator D and the parameter e(k) are given by the following formulas, respectively. 

D = e (fi(k)^ls^k,,+i^ljk,,), (2-9) 

-£-<fc<£:,s 



e(k) = y '[2t cos(ka)} 2 + [Asin(ka)} 2 , (2.10) 
The coefficient (3{k) in Eq. (2.9) is given by the formula: 

The coefficient /3(k) is related to the interband transition between the conduction and valence bands in a unit cell of 
length 2a, and the second term in Eq.(2.9) is often related to the intraband transition [13,14]. 



B. Berry phase and analytical format of one-dimensional Bloch functions 

In periodic systems, Bloch functions can not be analytical and are discontinuous for a composite band in which one 
band structure contains more than one branch. This fact was firstly pointed out by E.I. Blount [44], and later proved 
by J. Zak [42]. Generally speaking, Eq.(2.4) is not analytical for composite bands. Fortunately, for the ID periodic 
system being discussed here, it has been proved by Kohn that analytical results can be obtained [45] since both valence 
and conducting bands are simple. Thus, we may avoid the trouble due to the discontinuous or non-analytical feature 
of Bloch wave functions. 

The Berry phase for a ID crystal with the centro- or inversion symmetric only can only be or it (mod 2w) [46], 
therefore one can expects a vanishing Berry phase for closed path in this specific problem. The Berry-phase in the 
crystals is related to the diagonal matrix of first term in Eq.(2.9) [46,47]. However, in the problem we are discussing, 
the Berry phase is 0. 
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C. Field theory of hyperpolarizabilities 



There are a lot of different formulas to compute the hyperpolarizabilities - Orr-Ward sum-over-state (SOS) method 
[48], dipole formulas by Shen [49], and Genkin-Mednis approach [50] generalized by Bishop and his coauthors [29]. 
Different gauge approaches can be applied to this problem. All are based on the perturbation expansion. Recently, 
the field theory has been used to discuss this nonlinear problem [13-15]. 

For this single electron problem, the Feynman diagram of x (3) is simply described as one connected cycle in Fig.l. 
If the long-wavelength approximation is applied, ki = k 2 = k 3 = K = 0. The only permutation to be considered in 
this graph is three different frequencies w\, uj 2 and w 3 . 

The general third-order susceptibility under D • E gauge is described by: 



x (3) (^; wi,w 2 ,w 3 ) = 



l r - 3 



/ 



dri<ir2(Ir3 J dtidt 2 dt 3 
rfrdte- iK - r + im (fD(r,i)D(r 1 ,i 1 )D(r 2 ,t 2 )D(r3,t 3 )), (2.12) 
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where V is the total volume, £1 = — Y^Wj, T is the time-ordering operator, D is the dipole operator, and (• • • ) 

i=l 

represents the average over the unperturbed ground state. 

For periodic systems, in order to maintain the periodicity of the position operator r, a " saw-liked" position operator 
must be introduced [51]. For convenience in studying the nonlinear susceptibilities, we usually express Eqs.(2.12) in 
the momentum space for further calculations [13,14,43]. 



III. HYPERPOLARIZABILITIES FOR SSH AND TLM MODELS UNDER DIPOLE FORMULA 

The TLM model [32] is simply the continuum limit of the SSH model [31], and analytical results for the TLM 
model can be easily derived from those under the SSH model. Therefore in the following part, we will first focus on 
the hyperpolarizabilities under the SSH model. Then we will deduce results under TLM model by simply passing 
to the continuum limit. As for the notations \b, the superscript A represents the abbreviation for the different 
four-wave-mixing (FWM) terms, the subscript B represents the different models such as SSH or TLM models. We 
also use \ to represent the hyperpolarizabilities without considering \7 k term in Eq.(3.1). 

A. General four-wave-mixing(FWM) results 

Under SSH and TLM models, the general four-wave-mixing(FWM) can be expressed as XssH^i^ii w 3 )= 
Xssh(-(^i +u 2 +w 3 ),wi,w 2 ,w 3 ): 

„4 



d 

{P{k)cr 2 + i — )G(fc,w-wi) 
d 

(P(k)a 2 + i-^)G(k, w - wi - uj 2 ) 

(P(k)a 2 + i-j^)G(k,LJ - W1-U2- W3) j, (3.1) 

where L is the chain length, n is the number of chains per unit cross area, and V(u>i, w 2 , w 3 ) represents all permutations 
for u>i, u} 2 and uj 3 (therefore the intrinsic symmetry is maintained [39]). The polymer chains are assumed to be oriented, 
and Green's function G(k,w) is defined as follows [13,14,43]: 

uj 2 - LUt + Id 
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with ojk = e(k)/h and e = + . 

After tedious derivations in the Appendix B, we obtain the following analytical results for the SSH model: 

15 f 1/s xdx 



^/(1-S 2 x 2 )(x 2 - 1) 
(2x - z\ - z 3 ) 



(3.3) 



where 



and 



X 8 (x - Zx)(x + Z 2 )(X - Z 3 )(x -Z-L-Z2- z 3 ) 

(2x + z\ + z 3 ) 
x 8 (x + z\)(x — z 2 )(x + z 3 )(x + z\ + z 2 + z 3 ) 
4(1 - 8 2 x 2 )(x 2 - 1) (3x - 2zi)(3x - 2( Zl + z 2 + z 3 )) 

x 8 (x — z\ — z 2 ) (x - zi) 2 (x - z\ — z 2 — z 3 ) 2 
4(1 - S 2 x 2 )(x 2 - 1) (3x + 2 Zl ){3x + 2{ Zl + z 2 + z 3 )) 
x 8 (x + zi + z 2 ) (x + zi) 2 (x + Z\ + z 2 + z 3 ) 2 

(3) _ 8 e'n^af 
Xo ~ 45 irA^ [6A) 

Zi = fori = 1,2, 3. (3.5) 

By setting z\ = z 2 = z 3 = z = huj/2A, Eq.(3.3) can be simplified as third harmonic generation {— 3w; u), u), oj). It 
is easy to prove that Eq.(3.3) is the same as Eq.(2.20) [13] or Eq.(9) [14] in our previous works. 

By changing x — > x + ie in Eq.(3.3), and by choosing the same parameters used in our previous works for polyacety- 
lene [13,14,21,22], A = 0.9eV, n = 3.2 x 10 14 cm- 2 , a = 1.221 and e - 0.03, we have 5 = 0.18 and Xo 3) « 1.0 x lO" 10 
esu. The absolute value of FWM is plotted in Fig.l. 

From the graph, we find several symmetrical resonant frequencies. The biggest resonant peaks are around (0, 0, ±1), 
(±1, ±1, =Fl) an d their permutations, which correspond to DC Kerr effect(DCKerr) and Optical Kerr effect or 
Intensity-dependent index of refraction(IDIR), respectively. There are some secondary resonant frequencies, the cusps 
shown in the Fig.l are around (±1, =f1, 0) and their permutations. They correspond to DC-electric-field-induced op- 
tical rectification(EFIOR). The resonant peaks for third- harmonic generation (THG) are not obvious in this graph. 

As to the magnitudes of resonant peaks, we have x^ 3 H — w ! 0, 0, a;) » ( — w ; w , w , — w ) ^ X^(0j w, —to, 0) ^> 
X (3) (-3o>; u),u),u). 

Eq.(3.3) could be further simplified as hyperpolarizabilities under TLM model by letting 5 —* 0: 

X^(^ 3 ) = { ^(4, Z) - ^ M(4, Z) + ± ± £ ( 4, „) 

V- z 2 {-z l + 2(z 2 + z 3 )) 
+ P £ z) 2z 2 z 3 (z 2 + z 3 ) L ^ 

,^xi ] \ V iz l+ z 2 f( Zl+Z2 -2z 3 ) 
+ M(4,z 1 + z 2 ) 



(3.6) 



P(zi,Z 2 ,Z3) 



+ ^ lV ^ 1 ~^: 2 2 ;^" rU ^ 3; M(4,z 1 ) 



-zf (zf - 2zi(z 2 + z 3 ) + 6z 2 z 3 ) 

T 

where 



2^2 z 3 



a := (z 1 + z 2 )(z 2 + z 3 )(z 3 + zi), (3.7) 

Z := z x + z 2 + z 3 , (3.8) 
L(n, z) and M(n, z) (n = . . . 4) are defined by (A2) and (A3) respectively. 
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B. DC Kerr effect 



By setting z\ = z 2 = and z 3 = z (z = fiu/(2A)) in Eq.(B36), we have 

■ 1/s dx 



(_J2x 
\ x 2 (x 



^J{l-5 2 x 2 ){x 2 - 1) 
-z) (2x + z) 2 



x — z) 2 x 2 (x + z) 2 x(x 2 — z 2 ) 



+ 2{l-5 z x z ){x z -1) 



3(3ai-2z) 3(3x + 2z) 



(3.9) 



x 2 (x — z) 2 x 2 (x + z) 2 
(3x-2z) 2 (3a; + 2z) 2 3(3x - 2z) 3(3a; + 2z) 
[x — z) 5 (x + z) 5 x(x — z) 3 x(x + z) 3 

After simplifications, we obtain the DC Kerr effect(DCKcrr) coefficient under the SSH model: 

DCKerr ((] n x _ (3) 15 cfe 



2 f 



2 — z 2 ) 2 x s (x 2 - z 2 ) 



+ 2{l~5 2 x z ){x z -1) 



16 12 



(3.10) 



(x 2 — z 2 ) 5 x 2 (a; 2 — z 2 ) 4 



1 I 



x 4 (x 2 - z 2 ) 3 x 6 (x 2 -z 2 ) 2 x 8 (x 2 -z 2 ) 
As <5 — > in Eq.(3.10), we have the results under TLM model: 

1 r 

Xt C lm"\^ 0, u>) = X ( 3) ^ (-2L(3, z) - L(4, z) + 32M(0, z) n) 
+ 24M(1, z) + 2M(2, z) + 2M(3, z) - 6M(4, z)) , 
where L(n, z) and M(n, z) (n = . . . 4) arc defined by (A2) and (A3) respectively. 

Substituting the identities for L(n, z) and M(n, z) into the above equation and simplifying the resulting expression, 
we obtain 



(3.12) 



X?2m"^ 0,w) = X ( o } 512z8( * 5 2 _ 1)3 {-(120z 8 - 580z 6 + 1029z 4 - 780z 2 + 216)/(z) 

+ -^(384z 12 - 928z 10 + 760z 8 - 22182z 6 + 65541z 4 - 66780z 2 + 22680) 1 
105 J 

Dropping the terms related to Vfe or in Eq.(3.1), we have 

*Sgg' r *m<>>) - tf^l - wrSw^T) K 6 (- 2 -- 2 ) 2 " xtj-z 2 )} (3 - 13) 



As S — > in the Eq.(3.13), we have 



XrfM err (0, 0, u>) = J| (2L(3, z) + L(4, z)) , (3.14) 

where L(n,z) is defined by (A2). 

Substituting the identities for L(n, z) into the above equation and simplifying the resulting expression, we obtain 

1 , (3-15) 

+ Y?j^( 48z8 _ 104z6 _ 154z4 _ 315z2 + 63 °) r 
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C. DC induced second harmonic generation 



By setting z\ = z 2 = z (z = fou)/(2A)) and z 3 = in Eq.(B36), we have 

dx 



{ 



y/(l-6 2 X 2 )(x 2 -1) 

(2a; - z) (2a; + z) 



x(x 2 -z 2 )(a;-2z) a:(x 2 - z 2 )(a; + 2z) 
1 1 



+ an - s 2 t 2 Mt 2 - u i 3x( * 3x - Az) + 3x(3x + 4z) 

+ 1 ^ J1 z 2 (x-z)(a;-2z) 2 z 2 (x + z)(a; + 2z) 2 



(3.16) 



x(x — z)(x — 2z) x(x + z)(x + 2z) 

^-^-"{ ^(^ 

(3a:-2z)(3a:-4z) (3x + 2z)(3x + 4z) 
(a; - z) 3 (x - 2z) 2 (x + z) 3 (x + 2z) 2 
(3a; - 2z) (3a; - 4z) (3x + 2z) (3a; + 4z) 
+ (x - zf(x - 2zf + (x + z) 2 (x + 2z) 3 

After simplifications, we obtain the DC induced second harmonic generation(DCSHG) coefficient under SSH model: 

v BCSflG/ n . .\ 15 v (3) / 



128 A0 ^(l-^z 2 )^ 2 - 1) 

1 



+ 



x 8 (x 2 — Az 2 ) x 8 (x 2 — z 2 ) 
+ 2(l-^V)(x 2 -l)| « + 1 



(3.17) 



_x 8 (x 2 -z 2 ) x 6 (x 2 - z 2 ) 2 x i (x 2 - z 2 ) 3 
36 28 32 

~ .T 8 (x 2 _ 4z 2) + 3.6(3.2 _ 4z 2)2 + 3.4(3.2 _ 4z 2)3 

As 5 — > in Eq.(3.17), we have the result under the TLM model: 

Xrf^ G (0, w, w) = ^^ 3) (-4L(4, 2z) + L(4, z) + 12M(4, z) + 2M(3, z) ^ 
-8M(2, z) - 72M(4, 2z) + 56M(3, 2z) + 64M(2, 2z)) , 

where L(n, z) and M(n, z) (n = . . . 4) are defined by (A2) and (A3) respectively. 

Substituting the identities for L(n, z) and M(n, z) into the above equation and simplifying the resulting expression, 
we obtain 

X^ G (0, u, u) = xP 256z8 ™ 2 _ 1)2 {(28z 6 - 104z 4 + 118z 2 - 43)/(z) 



- tL( 288210 - 184 ^ 8 + 204 ^ 6 - 5068z 4 + 9380z 2 - 4515) 1 
105 J 

+ ^ 3) 2048z 8 ?4 5 z 2 -l) ^- 48 - 4 + fe2 + 1 ^ 



(3.19) 



+ — (12288z 8 + 4992z 6 + 2464z 4 - 910z 2 - 105) 1 
105 J 

Dropping the terms related to or ^ in Eq.(3.1), we have 
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As d -> in Eq.(3.20), we have 



x£££ HG (0, W ) = X < 3) (-4L(4, 2z) + L(4, z)) , (3.21) 
where L(n,z) is defined by (A2). 

Substituting the identities for L(n, z) into the above equation and simplifying the resulting expression, we obtain 

X^ G (0,c, co) = X^Y^g {s.f (*) - \f(2z) -i(128z 4 + 200z 2 + 315)} (3.22) 

D. Optical Kerr effect or intensity-dependent index of refraction 

By setting Z\ — z 3 = z (z = 7kj/(2A)) and z 2 = — z in Eq.(B36), we have 

dx 



I" 



a;V(! - <5 2 z 2 )(z 2 - 1) 
4x 1 1 



(x 2 -z 2 ) 2 (x-z) 3 (x + z) 3 

(3x-2z) 2 (3x + 2z) 2 



(3.23) 



+ 2(l- ( 5 2 x 2 )(x 2 -l) 1 , 

x(x — zy x(x + zp 

(3x-2z) 2 (3x + 2z) 2 2(9x 2 -4z 2 ) 

+ (x - 2z)(x - z) 4 + (x + 2z)(x + z) 4 + x(x 2 - z 2 ) 2 

After simplifications, we obtain the optical Kerr effect or intensity-dependent index of refraction (IDIR) coefficient 
under the SSH model: 

idir, \ (3) 15 f 1/S dx 

Xssh (w.-w.w) =Xo 



128 ^ ^(l-^x 2 )^ 2 - 1) 
1 4 



x 6 (x 2 -z 2 ) 2 x 4 (x 2 -z 2 ) 3 
+ 2(l-<5 2 x 2 )(x 2 -l)(- 10 18 



(3.24) 



x 8 (x 2 - z 2 ) x 6 (x 2 — z 2 ) 2 
64 



x 4 (x 2 — z 2 ) 3 x 8 (x 2 — 4z 2 ) 
As 5 — > in Eq. (3.24), we have the result under the TLM model: 

X^if (w, -w, w) - x[, 3) ^ (i(3, *) - 4L(2, z) - 32M(4, z) ^ 
-26M(3, z) - 16M(2, z) + 128M(4, 2z)) , 
where L(n, z) and M(n, z) (n = . . . 4) arc defined by (A2) and (A3) respectively. 

Substituting the identities for L(n, z) and M(n, z) into the above equation and simplifying the resulting expression, 
we obtain 



IDIR (3) 15 f 8z 8 + 44z 6 + 83z 4 - 495z 2 + 315 

tlm -X 25^| 15(z 2 -l) 2 



4z 6 + 22z 4 - 43z 2 + 20 ,, , 

-/(z) + (4z 2 - l)f(2z) 



(3.26) 



(z 2 - l) 2 

Dropping the terms related to Vfe or in Eq.(3.1), we have 

&% R (u>,-u,u) = A 3) — f 1/S - { - - 1 (3.27) 

XssH[ > ' > Xo 128 A 7(1— I^X^T) \x6(x 2 -z 2 ) 2 x 4 (x 2 -z 2 ) 3 / ^ U 
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As 5 — > in the Eq.(3.27), we have 

x'tTm^, = (i(3, *) - 4L(2, z)) , (3.28) 

where L(n,z) is defined by (A2). 

Substituting the identities for L(n, z) into the above equation and simplifying the resulting expression, we obtain 

X™'«(0, u,u) = ^ 3) 25fe8(z 1 2 _ 1)2 {-(840z 4 - W25z 2 + 630)/(z) 
+ (16z 8 + 88z 6 + 226z 4 - 1005z 2 + 630)} 



E. DC-electric-field-induced optical rectification 



By setting z\ — 0, z 2 = —2, z 3 = z (z — fiw/(2A)) in Eq.(B36), we obtain the DC-electric-field-induced optical 
rectification (EFIOR) under the SSH model: 

15 f 1/s dx 



EFIOR/r, \ (3) 15 /" 

Xssff lU, -w, — Xq ^ J 



X 



5 2 x 2 )(x 2 -1) 



3x 2 + z 2 2 2 2 2z 4 + 9x 4 - 7x 2 z 2 

+ 6(1 — o x )(x - 1) 



v (»)il /" _ (3.30) 
Xo 128 A 



(x 2 — z 2 ) 2 ' (x 2 — z 2 ) 3 

^ dx_ 

^(1-S 2 x 2 ){x 2 -1) 
-2 1 
x 6 (x 2 — z 2 ) 2 x 8 (x 2 — z 2 ) 

+ 6(1 - S 2 x 2 )(x 2 1) f + 3—2^—^2 + 



x 4 (x 2 -z 2 ) 3 x 6 (x 2 -z 2 ) 2 x 8 (x 2 -z 2 ) 



As 6 — > 0, we have the result under the TLM model: 



xf -w, w) = X^ 3) (-2L(3, z) - L(4, z) + 6(4M(2, z) + 3M(3, z) + 2M(4, z))) , (3.31) 
where L(n, z) and M(n, z) (n = . . . 4) arc defined by (A2) and (A3) respectively. 

Substituting the identities for L(n, z) and M(n, z) into the above equation and simplifying the resulting expression, 
we obtain 

X§El? R {0-u>,v) = X^ 3) ^Li) {-W5 {2AzS 38z6 + U9Z * ~ 252 ° z2 + 2520) (3.32) 
+ (15z 4 -40z 2 + 24)/(z)} 
Dropping the terms related to Vfc or ^ in Eq.(3.1), we have 

^°*( '- w ' w >=xS 3) ]i/ ^(1 - g^x* 2 - 1) i * 6 (* 2 2 ^ 2 ) 2 " z2 ) 1 (3 ' 33) 

As 5 -> in the Eq.(3.33), we have 

Xt lm *(0, -w, w) = x[, 3) ("2£(3, *) - i(4, z)) , (3.34) 

where L(n, z) is defined by (A2). 

Substituting the identities for L(n, z) and M(n, z) into the above equation and simplifying the resulting expression, 
we obtain 

y$l$ R %-w) = xi 3 \ 28 J 5 a _ n {(7* 2 - 6)/(z) 

1 , (3-35) 

+ y^( 48z8 _ lMz& ~ 154z4 _ 315z2 + 63 °) 
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F. Third harmonic generation 



The results of third harmonic generation(THG) can be computed by setting z\ = z 2 = 23 = z in Eq.(B26). They 
can also be found in our previous works [13,14]. Here we simply state the result: 



JTHG, x (3) 45 H 

Xssh i")-Xo Y28 A 77r= 



dx I 47 - 48(1 + (5 2 )x 2 + 48c5 2 x 4 



vtnv)Fi)l 8x8(x 2 -z 2 ) 

3(l-<5 2 x 2 )(x 2 -l) 

x 6 (x 2 — z 2 ) 2 
9 [47- 48(l + (5 2 )a 



+ 
+ 



8x 8 (x 2 - (3z) 2 ) 

63 (l-£ 2 x 2 )(x 2 - 1) ' 
x 6 (x 2 - (3z) 2 ) 2 



and 



Excluding Vfc terms, we obtain 

X SSff Xo / v/(1 _ (52x2)(x2 _ 1) ( ,8 (x2 _ z2) - ^ _ (3z)2) ) (3-38) 
As S -> in the Eq.(3.33), we have 

w, w) = X^g^I {5(729/(z) - /(3s)) - 8(216z 4 + 300z 2 + 455)} (3.39) 



IV. DISCUSSIONS 



A. resonant and non-resonant features 



Fig. 3 shows the hyperpolarizabilities of DCKerr, DCSHG, IDIR and EFIOR for SSH models. Choosing the same 
parameters in Section III A, we see the magnitudes of the resonant peaks are in the following order: 

X™(-u;0,0,u>) > xS ) s (-^;^-^^)>xfe ) s (-2w;0,w,a;) 

> xill(0;u;-u,0)> xHI(-3u;lu,lu,lu). 

Fig.4 shows the comparison of hyperpolarizabilities of DCKerr, DCSHG, IDIR, EFIOR, THG for SSH models for 
low frequencies. Choosing the same parameters in Section III, we see the non-resonant features are in the following 
order: 

/o\ /o\ (r>\ 

XrJ™-re S (-3w;W,W,w) > Xnon-res (~ 2w ; 0, U), Uj) > Xnon-res (~^: ^ ■, ~U),Ul) 

(3) (3) V^-^> 

> X non—res 

(-w;0,0,w) > x 

non — res 

(0; ui, -w,0). 

This non-resonant frequency dependence of nonlinear optical properties is consistent with previous calculations for 
polyenes [27]. 
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B. Vfc term, Krammer-Kronig(KK) relation, the overall permutation and Kleinman symmetries 



There are some arguments about whether or not one should include the terms in the calculations of nonlinear 
optical properties. The terms are usually considered to be related to the intra-band current [12]. Otto suggested 
not to include this term because of the non-periodic property of V& [28,34,36], and thus the calculations would be 
purely based on the inter-band transition. In this work, we compute the analytical results with or without Vfe terms 
to show the differences of the results under both schemes. 

Considering the important physical contributions from the intra-band currents [12,29], we are in favor of including 
the Vfe terms. Moreover, the restriction of our calculations in a unit cell implicitly imposes the periodic condition 
even for Vk operator. For the linear susceptibility x^\ our calculations [43] show that the Vfc term make no 
actual contributions. Due to the centro- or inversion symmetry for both SSH and TLM models, the second-order 
susceptibility x^ vanishes, the first nonzero susceptibility is the third-order susceptibility x^ ■ From the formula, 
the Vfe term causes the non-commuting problem between operators in the third-order susceptibility x^ calculations. 
Therefore, it breaks the overall permutation symmetry (between (wi, u>2, 0J3) and f2 = — that is preserved in 

molecular systems [39] where only bound states exist. Kleinman symmetry, which is defined as the intcrchangcability 
of all n indices in the rank n tensor x^ [37,39] and derived from the overall permutation symmetry, is also broken 
even for low optical frequencies. Our calculations clearly show the nonequivalcnt off-resonant behavior between 
X DCKerr (—uj; 0, 0, u) and x EFIOR {0] — w, 0). To provide a theoretical result that can be measured by experiments, 
we also perform the calculation of x^ ( w ; w i w > — 3w), which is the overall permutation of the THG x^ ( — 3w; to, w, u). 
Wc obtain: 



{ 



l ' s dx 



y/(l ~5 2 X 2 ){x 2 - 1) 

27 



x 8 (x 2 - z 2 ) x 8 (x 2 -9z 2 ) 



+ 8(l-6 2 x 2 ){x 2 -1) 



74 1 11 1 



(4.3) 



9x 8 {x 2 -z 2 ) 3x e {x 2 -z 2 ) 



2 



512 1 1 1 U 

+ ~9~x s (x 2 -4z 2 ) ~ M x s {x 2 -9z 2 ) + 63 x 6 (x 2 -9z 2 ) 2 J J 

Letting S — > in (4.3), we obtain the result under the TLM model: 

xP LM (co;co,co,-3co) = X^^E {^z 2 61)/(z) + ^(4z 2 - l)/(2z) 

"243 (1944z2 " 241 ) / ( 3z ) + 243 (27z4 - 30z2 + 805) | (44) 

13) , 5 80 , 28500 , ^, , 
= ^28 + 33* + TooT* +0{Z)) ( ^ 0) ' 

On the other hand, analytical expression for THG under the TLM model is given by the formula (see, for example, 
[13,14]): 



Y (3) ( o w . w W ^- Y (3) 45 i 14 4 (37-24z 2 ) (lz**l mz) 
Xtlm(-^"^) - Xo 128 1 -3^8 " 15? + §P /( } + 24z« /(3z) 




(4.5) 



z A + 0(z b )) (s->0). 



("i) (3) 

Both x^| M (w;a;,w,— 3u) and Xtlm(~ 3w; uj, w, ui) in the non-resonant region are plotted in Fig. 5. The numerical 
calculation shows 40% difference between them when z = 1/6 (corresponding to 0.3eV or wavelength 4.14/zm) in 
this model. Thus it clearly shows the break of overall permutation symmetry, while the Kramers-Kronig relation 
is satisfied in Eq.(4.4) and Eq.(4.5). Subtracting two asymptotic expressions in Eq.(4.4) and Eq.(4.5), we have the 
following relationship of the difference in the off-resonant region: 

(3) e 4 n 4a 3 h 2 uJ 2 
5 X (i, (uj)ot — . (4.6) 
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By excluding the Vfc terms, our calculations show both the overall permutation and Kleinman symmetries remain 
valid. For example, x DCKerr =x EFIOR an d x^ (— 3^; w, w, w)= w, w, — 3w) are preserved for all frequencies. 

Obviously, it is the \7 k term that breaks the overall permutation in periodic systems. Though the experimental 
stimulated Raman scattering does not exhibit the overall permutation symmetry for the resonant region [38,39], it 
seems that no experiment has been reported to test the validity of the overall permutation symmetry for low frequency 
and non-resonant region. Therefore, we suggest a new x' 3 ^ experiment on infinite trans-polyacetylene chains to test 
the break of overall permutation symmetry. 

Recent experimental studies have already pointed out the deviation from the Kleinman symmetry for low optical 
frequencies in numerous investigations of various optical systems [38] . The assertion of the general failure in Kleinman 
symmetry has been made in [38]. Due to the similarity between the overall permutation and Kleinman symmetry, 
Eq.(4.6) can also be used to explain the break of Kleinman symmetry in experiments qualitatively. Eq.(4.6) shows that: 
(i). The break increases with decreasing band gap and is proportional to J 1 . These are consistent with the previously 
reported experiment [52]; (ii). The break increases with to (the hopping of 7r electrons between the nearest- neighbor 
atoms). This explains the experimental results that the deviation of Kleinman symmetry is favorable of (20^50%) the 
delocalized states such as aromatic molecules [53] and some polymers or crystals [54,55], while unfavorable of(<8%) 
the localized states such as molecular systems such as 02, N2, etc. [56,57]. On the other hand, the vanishing x^ 
under the SSH or TLM model shows that some symmetries such as centro-symmetry, etc, can suppress the deviation 
from Kleinman symmetry even for periodic systems. This may explain why Kleinman symmetry is still preserved in 
some x^ experiments of crystals [58]. 

The magnitudes of the hyperpolarizabilities with terms are quite close to those without Vfc terms in our results. 
In this sense, both results can give the correct position of the resonant peaks qualitatively. However, we should 
notice that there is actually a sign difference between the results of including and excluding Vfc terms. Therefore, the 
intraband current contribution to the hyperpolarizabilities is actually about as twice big as the interband contribution 
and can not be neglected. The total contribution (without terms), interband contribution (with Vfc) and intraband 
contribution (the difference of previous two) under TLM model are plotted in Fig. 6. 



C. Zero frequency behaviors for hyperpolarizabilities 

Hyperpolarizabilities under zero frequency arc also called static hyperpolarizabilities. To obtain the static hyperpo- 
larizabilities, there are several different approaches. One way is to obtain the static polarizability by directly applying 
the static electronic field. For example, by applying Rcsta's definition of dipolc moment for periodic systems [59], 
Soos and his coauthors obtained the polarizability of one-dimensional Peierls-Hubbard model based on the static 
electronic field [60]. The other way is to obtain the optical hyperpolarizabilities first, then let all frequencies approach 
0. Obviously, both ways should yield the exactly same results of static hyperpolarizabilities. Here we applies the 
latter method to do our calculations. 

We may use the results in Section III under TLM model to study the zero-frequency behaviors, with or without 
Vfc terms. The zero-frequency results under SSH models shall be similar. 

1. DCKerr 

Let z — > in Eq.(3.12), we obtain 

XS-(0,0,- 0) = x^ 3) (| + + + ^ + 0(z*)) (4.7) 



Dropping Vfc term and let z — > in Eq.(3.15), we obtain 

-DCKerr / n n ,, _^ n \ (3), 1 50 2 40 4 48 g g 

Xtlm (0,0,w^0)=x ~ Y43 ~ 143 >> t ' 



2. DCSHG 

Let z — > in Eq.(3.19), we obtain 

ocshg (() u u n] <3) ( 5 254 68620 4 375824 
Xtlm (0,w - 0,u - 0) - Xo (^ + ^f z + ^oT + ^003~ Z + ° {z )} (49) 

Dropping Vfc term and let z — > in Eq.(3.22), we obtain 

X^ G (0,^ 0,c^ 0) = X f{-\ - f f 2 ^/ - ^ + 0(z*)). (4.10) 
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3. IDIR 

Let z — > in Eq.(3.26), we obtain 

xKifta, - 0, -a, - 0, w - 0) = X ^ 3) (| + |z 2 + + ^z 6 + 0(z 8 )). (4.11) 

Dropping Vfe term and let z — > in Eq.(3.29), we obtain 

xSi?( W - 0, -a, - 0,« - 0) = x^(-7 - ^ ^ ^/ + 0(z 8 )). (4.12) 

4. EFIOR 

Let z — > in Eq.(3.32), we obtain 

xfflT(0, - w - 0, w - 0) = x 3) (| + g* 2 + ^z 4 + + 0(z*)). (4.13) 



Dropping Vfe term and let z — > in Eq.(3.35), we obtain 

xS°>, -u,,u,) = X f \-\ ^ - ^z 4 - ±Lf + 0(z 8 )). (4.14) 

5. THG 

From Eq.(3.3) in our previous work [13], let z — > 0, we have 

T „ r/ „ „ „. (3), 5 80 n 98580 , 96640 fi ^. 8 ., ,„ , 

X?f^(u, - 0, w - 0, u -> 0) = xf, 3) ( + rf z + + + 0(z 8 )). (4.15) 

Dropping Vfe term and let z — > in Eq.(3.39), we obtain 

*??£(a, - 0,. - 0,0, - 0) = x< 3) (-i - ^z 2 - ^ i|f. 6 + 0(. 8 )). (4.16) 

From the zero frequency behaviors, we can see the difference of results between having and not having V& terms. 
Although the magnitude of hypcrpolarizabilities at (0,0,0) is quite close, 5x /28 vs. 4x /28. 



D. Comparison with other theoretical results 

For the nonlinear properties under single electron periodic models, Genkin-Mednis developed an approach [50] that 
was later applied to polymer systems by Agrawal, et al [12]. In those works, general formulae of nonlinear optical 
response were developed. BY applying the Genkin-Mednis approach to the SSH model, Wu and Sun obtained the 
analytical format of third-harmonic generation (THG) [21,22] that yields the same results as those under the static 
current-current ( Jo Jo) correlation [15]. As we discussed previously [13,14], the general results obtained here are quite 
similar but different from the results obtained before. On the THG problem, our results are qualitatively close to 
the numerical results obtained by Yu, et al [19,20] and Shuai, et al [23] 's works. The reason for those difference is 
caused technically by the treatments of Vfe operator, but physically by the gauge phase factor [43] that we are going 
to discuss further in the subsequent paper [40]. 

Recently Kirtman, Gu and Bishop extended Genkin-Mednis method to the fully coupled perturbed Hartree-Fock 
theory in discussing the hyperpolarizabilities of infinite chains [29]. In our results, it is still an uncoupled treatment 
because the fully-filled valence and empty conduction band structures in both SSH and TLM models lead no difference 
in the final results [29] . 

Finally, we would like to point out that the results under DD correlation are also different from those under J Jo 
correlation [15,16]. More details will be discussed in our subsequent paper [40]. 



E. Some suggested experiments 

Our calculations predict the break of overall permutation symmetry of hyperpolarizabilities for periodic systems in 
off- resonant regions. Despite of the wide acceptance of overall permutation symmetry in off- resonant regions [39], the 
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direct measurement of the above assertion does not appear to be present in the literature. In our previous work [61], 
we suggested experimentalists perform off-resonant measurements on some centro-symmetric ID periodic systems 
such as trans-polyacetylene, etc to directly test the break of overall permutation symmetry and Eq.(4.6). Though 
our calculations are based on ID periodic models, the results may be applied to some delocalized two-dimensional 
materials such as benzene ring structures where the periodicity is along the circle of ring. Therefore, experimentalists 
can also perform some off-resonant x^ experiments on some symmetric aromatic molecules to directly test the overall 
permutation symmetry of hypcrpolarizabilities. 

Besides the break of overall permutation symmetry, the relationships between different nonlinear optical processes 
such as Eq.(4.1) and Eq.(4.2) can also be tested experimentally. 



V. CONCLUSIONS 



Analytical expressions of the general four-wave-mixing (FWM) for infinite chains under the SSH and TLM mod- 
els are first derived under the DD correlation by the field theory. The hyperpolarizabilities of DCKcrr, DCSHG, 
IDIR and EFIOR are obtained by either including or excluding terms. Contrary to the linear response, in- 
traband contributions to the hyperpolarizabilities dominate the final results if they are included. Vfe term also 
breaks the overall permutation and Kleinman symmetry in even for the low frequency and non-resonant re- 
gions. The reported deviations from Kleinman symmetry in experiments are explained qualitatively in this work. 
A new feasible off-resonant x*- 3 ^ experiment is suggested to test the break of overall permutation symmetry in in- 
finite ID periodic polymer chains with centro-symmetry such as trans-polyacetylene, etc. For the infinite single 
electron periodic systems, our calculations show the following trends for the various third-order nonlinear optical 

/o\ /o\ /o\ 

processes in the non-resonant region: Xnon—resi — 3cj] cj : cj : cj) > Xnon—res 

XnL-res(-u; 0, 0, w) > XnL-res( ; w > °)> and in thc resonant region: xr? s (- w ;0,0,w) > Xr%(-u;uj,-(j,u) > 
Xres(— 2w; 0,lu,lu) > Xrel (0; uo] — oj , 0) > Xrel (— 3w; u), u), lo). We look forward to the experimental testings on thc above 
theoretical results. 

This analytical calculations of hyperpolarizabilities are tedious, but helpful to illustrate both single electron models 
and theoretical methodologies in nonlinear calculations. Based on single electronic models, more sophisticated models 
like Hubbard model [17,24,26], electron-hole pair model [17,25], etc, may be studied by applying similar techniques 
developed here. 



APPENDIX A: SOME INTEGRALS 



Define 



/(*) := 



dx 



i (x 2 - z 2 )Vx^l 



arcsin(z) 
cosh _1 (^ 



zy/z 



2zVz 



(z 2 < 1), 
(z 2 > 1). 



(Al) 



r 

L(n,z) := J 



dx 



r 2n 



(x 2 - z 2 ) 5 -«v / ? r— T 



(A2) 



\/x 2 — ldx 
x 2n (x 2 - z 2 ) 5 ~ n ' 



Using Maple, we obtain 



L(4, z) = (105/ (z) - (48z 6 + 56z 4 + 70z 2 + 105)) , (A4) 
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i(3, z) = 30^8^2,^ ((-120z 2 + 105)/(z) + (16z 6 + 24z 4 + 50z 2 - 105)) , 



L ( 2 < *) = 24,8(^2 _ 1)2 (( 144 ^ 4 - 240z2 + 105 )/W 
- (16z 6 + 40z 4 - 170z 2 + 105)) , 

M(4, z) = ((105z 2 - 105)/(z) - (8z 6 + 14z 4 + 35z 2 - 105)) 



M(3, z) = ^ ((-90z 2 + 105)/(z) + (4z 4 + 20z 2 - 105)) 



M (2, z) = 24z8(z2 _ ^ ((72z 4 - 180z 2 + 105)/(z) + (-8z 4 + 110z 2 - 105)) 



1 

z^ 

+ (-92z 4 + 200z 2 -105)) 

M( °' Z) = 384z»(^ - l) 3 ((_192Z6 + 432Z ' " 36 ° Z2 + 105)/(Z) 
+ (48z 6 - 248z 4 + 290z 2 - 105)) 

APPENDIX B: THE DERIVATION OF FWM UNDER SSH MODEL 

We're trying to calculate the following integral: 

X(u u U2, w 3 ) = - J 7^Tr(D k G{Q)b k G{LQi)i) k G(LUi + ^ 2 )^ fc G( Wl + lu 2 + w 3 )), 

where 

d 

D k := /3(fe)<7 2 + i— 

and 

w — s + Wj,cr 3 



(w - a) 2 —LJ k + ie' 
Let's define 

P± : = ^( 1±cr 3), 
i fc := /3(fc)cr 2 , 

4 := 4' 
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Then it's easy to prove that 



9+(s) ■■= 
9-(s) ■■- 



1 



uj — s — LUk + ie 
1 

LO — 8 + LO k — ^ 



G(s) = 9+(s)p+ + g-(s)p- 
A k p± =p T A k , 
BkP± =p±B k . 



1. 



J(wi, uj 2 , w 3 ) = h + h + h + h + h + h + h + h 



where 



2. 



dui 



TTr(ifcG(0)i fc G(wi)i fc G(wi + w 2 )i fc G(wi + lj 2 + lj 3 )) 



2^i 
duo 



: Tr(A fc G(0)A fc G(wi) J B fc G(wi + w 2 )B fc G(wi + lu 2 + lo 3 )) 



— Tr(A fe G(0)S fe G(a; 1 )^ fe G(a; 1 + w 2 ).B fe G( Wl + io 2 + w 3 )) 
2m 

- J ^-Tr(A k G(0)B k G(uJi)B k G(uji + w 2 )A fe G( Wl + uj 2 + w 3 )) 



:Tr(B fc G(0)i fc G(wi)i fe G(wi + W2)BfcG(wi + w 2 + w 3 )) 



dw 
2^i 

^ 7 Tr( J B fc G(0)i fe G(o;i)S fe G(a;i + w 2 )i fe G(^i + w 2 + w 3 )) 

27TI 



2^i 



Tr(B fc G(0)B fc G(wi)i fc G(wi + w 2 )i fc G(wi + w 2 + w 3 )) 



- y ^Tr(B fc G(0)B fc G(wi)B fc G(a;i + w 2 )B k G(u x +to 2 + w 3 )) 



^Tr(i fe G(0)ifeG(a;i)i fe G(a;i + w 2 )i fe G(u;i + w 2 + w 3 )) 
^— 5+(0).9-(wi).g+(wi + w 2 )p-(wi + w 2 + w 3 ) 
^— 5-(0)ff+(^i)5-( w i +^2).9+(^i +w 2 +w 3 ) 

= /l(Wfc) - Ji(-Wfc), 



h{u k ) 



f3 i {k){Auj k -uj 1 -uj 3 ) 



(2uj k - wi)(2uh + oj 2 )(2uj k - uj 3 )(2uj k -U1-U2- uj 3 ) ' 



h = ~ I ^Tr{A k G{Q)A k G{uj 1 )B k G{uj 1 +uj 2 )B k G(u 1 +u J2 +uj 3 )) 
J 2m 

/' duj d d 

— /3(%+(0)/3(%_(wi)— g-{tox +co 2 )—g-(u! 1 +oj 2 +w 3 ) 

+ J ^/3(%-(0)/3(%+(wi)^5 + (o;i+W2)^5+(wi+W2+w 3 ) 
= I 2 {to k ) - I 2 (-(J k ), 



(B2) 



(B3) 



(B4) 



(B5) 
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where 



3. 



where 



4. 



where 



d 



1 



)\uj=u> k ■ 



2u) k — U>1 dk LU — LUl — LU2 + UJk 9k LU — LUl — Ll>2 — Ll>3 + UJk 

/duj 
—Tr{A k G{<S)B k G{^)A k G{^ + W2)BfcG(wi + ^2 + W3)) 

= y ^r/?(%+(0)^:ff+(wi)/3(fc)g_(wi + w 2 )^5-(wi +^2+^3) 

+ / ^i^ 9 -^^ 9 '^ 1 ^^ 9 ^ 1 + ^§k 9+<yUJl +UJ2 + ^ 
= h(uk) - h(-Uk), 



h(uj k ) = 0{k) 



d 



0(k) 



d 



1 



<9fc (w — 0J2 + w/j)(a; — — u>2 + Wfe) dk uj — uji — 0J2 — UJ3 + uJk ^ Uh 

d (3{k) d 1 

dk u — 1V2 + u k dk (u) — u>2 — ^3 + Wfe)(o> — uj\ — U02 — ^3 + u k ) u ~ ulk 

dm k (3 2 {k) 1 d 1 

dk 2(2u> k — L02) 2u> k — L02 dk (2ujk — u>2 — u>3){2uj k — u>i — u>2 — L03) 

duj k P 2 (k) 1 d_ 1 

dk 2{2u k — L02) (2ojk — uj 2 )(2ujk — lu\ — lu 2 ) dk 2cj k — wi — u>2 — L03 

dm k (3 2 (k) 1 d 1 

dk 2(2u k — L02) 2ui k — U02 - ^3 dk (2uu k — co 2 — ^3)(^k - u>i — u) 2 - W3) 



du; k I3 2 (k) 



d 



dk 2(2uj k — U02) (2uj k - UJ2 — uj 3 )(2u; k — ui - 002 — 0J3) dk 2uj k — U02 - ^3 

-^rTr(A k G(0)B k G(LJi)B k G(u>i + w 2 )ifcG(u;i + lo 2 + w 3 )) 

= / ^i^ 9+ ^^ 9+ ^d~k 9+<yUJl + w 2)/?(%-( w i + ^2 + W3) 

/duj d c) 

^— 0(k)g-(O)—g-(uj-i)—g-(uJi + w 2 )/3(fc)s+(wi +uj 2 +^3) 

= h(<jJk) - h(-Uk), 



h(ui k ) = -P(k) 



d 2 



dk 2 (2uj k -W1—U2 — W3)(2uj k —UJ2 — oj3)(2uj k - 0J3) 
(3(k) d (3{k) d 1 



+ 



2 dk 2ui k - 


- L03 dk (2uj k — u>i 


- U>2 - U>3)(2u> k -U>2~ 


- w 3 ) 


P(k) d 


(3(k) 


d 1 




2 dk {2u k 


- u) 3 )(2uj k - lu 2 - 


L03) dk 2uj k — u>i — u>2 


- L03 


p 2 (k) d 


1 


d 1 


— )\ 



(B6) 



(B7) 



(B8) 



(B9) 



(BIO) 



UJ = UJ k ■ 
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5. 



where 



G. 



where 



7. 



where 



duj 



2m 
dui d 



T Tr(B fe G(0)i fe G(w 1 )i fe G(u; 1 + u 2 )B k G(uj 1 + lo 2 + w 3 )) 



d 



2ni (0)/?(%+ (wi )/?(%- (wi + u 2 )— + u 2 + w 3 ) 
+ y ^^5+(0)/3(%-(wi)/3(%+(^i +w 2 )^g+(cji + w 2 + lu 3 ) 
= h(u)k) - h{-Uk), 



d 



P 2 {k) 



d 



dk 2(2uj k + tJi)(2wk — w 2 ) dk 2ujk — oj 2 — oj 3 



/du> " 
-^—Tr(B k G(0)A k G(LJi)B k G(u;i + lu^A^^ + uj 2 + ws)) 

= y ^-^9-{0)P(k)g + (uJi)-^g + (uJi + uj 2 )(i(k)g-(uJi + lu 2 + uj 3 ) 

+ J ^^5+(0)/3(%-(wi)^fir_(a;i+a; 2 )/3(% + (a;i+a; 2 +a;3) 



i) 



d 0(k) 



dk (2uj k + u>i)(2u>h + oji + lu 2 ) dk 2u k - uo 3 

d p(k) d p(k) 

dk 2u> k + uji dk (2uj k — uj 2 — uj 3 )(2oj k — oj 3 ) 

d_ (3 2 (k) d 1 

dk 2(2uj k + L0i)(2uj k - lu 3 ) dk 2uj k - uj 2 - u> 3 

d_ (3 2 (k) d 1 

dk 2(2bj k + LUi)(2uj k — uo 3 ) dk 2ui k + u)\ + ui 2 



h = - / ^Tr{B k G(Q)B k G{Lu 1 )A k G(u 1 +u J2 )A k G(u Jl +u J2 +u J3 )) 



I 



duo d d 

2^idk 9 ^^dk 9 ^^ LJl ^^ 9+( '- LJl + w 2)/?(%- (^i + w 2 + w 3 ) 
+ ' ^^ff+(°)7^;5+(wi)/3(%-(wi +w 2 )/3(fc).g+(wi + w 2 +w 3 ) 
= I 7 (uo k ) - I 7 (-uo k ), 



I 7 (uo k ) = - 



d 



p\k) 



d 



1 



dk 2(2uo k + uo 2 )(2uj k — uo 3 ) dk 2uj k + uj\ + uj 2 
d 1 d (3 2 (k) 

dk 2uj k + ujx + uj 2 dk (2uo k + uo 2 )(2uj k - ui 3 ) ' 



(Bll) 



(B12) 



(B13) 



(B14) 



(B15) 



(B16) 
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8. 



h = ~ j ^Tr(B fc G(0)S fc G(o;i)B fc G(a;i +W2)BfcG(wi +wa +w 3 )) 

/d(jO 
—B k g + (0)B k g+(uJi)Bkg+(uJi + w 2 )Bfe.g+(wi + w 2 + w 3 ) 

—B k g-{Q)Bkg-{wi)B k g-{ui + w 2 ).B fe .9_(a>i + w 2 +w 3 ) 

= 

The expression for X{w\, u> 2 , w 3 ) is extremely messy. However, the quantity which we're really interested in is the 
sum of its six permutations 

S = ^C^ 1 '^'^) + Z(w 1 ,w 3 ,w 2 ) + J(w 2 ,wi,a; 3 ) 

6 (B18) 

+ I(lu 2 , w 3 , wi) + J(w 3 ,wi,w 2 ) +X(w 3 , w 2 , wi)) 

Due to its highly symmetric character, we expect a great simplification would occur. This is indeed the case as 
shown below. The key technique used is that we can permute u>i, oj 2 and w 3 freely when computing the summation. 
We present some details of the simplification procedure in the following. 



1. 



T , * T , v T , n d /3 2 {k) d 

h(u k ) + h(u) k ) + I 7 (uj k ) = - 



+ 



dk 2(2uj k + wi)(2wfe - lu 2 ) dk 2uj k - uj 2 - w 3 
d P{k) d f3{k) 



+ 



= + 



dk (2u> k + ui)(2uj k + loi + uj 2 ) dk 2uj k — lo 3 

o m o m 

dk 2oj k + u>i dk (2ui k — lo 2 — uj 3 ){2uj k — uj 3 ) 

d_ f3 2 {k) d 1 

dk 2(2uj k + u>i)(2uj k — L03) dk 2u> k — uj 2 — cj 3 

d_ f3 2 {k) d 1 

dk 2(2uj k + u\){2uo k — w 3 ) dk 2uj k + oj\ + uo 2 

d_ (3 2 (k) d 1 

dk 2(2uj k + u 2 )(2uj k — w 3 ) dk 2uj k + u)\ + uo 2 
d 1 d I3 2 {k) 

dk 2uj k + ujx + u 2 dk (2u k + uj 2 )(2u k - uj 3 ) 
d (3 2 (k) d 1 

dk (2uj k + ui){2uo k — uo 2 ) dk 2uj k — u 2 — uu 3 

d_ p(k) d (3(k) 

dk (2u k + LJi)(2u> k + loi + oj 2 ) dk 2uj k - lo 3 

d (3{k) d (3{k) 

dk 2oj k + uji dk (2co k — lo 2 — Lo 3 )(2u> k — w 3 ) 
d 1 d /3 2 {k) 

dk 2uj k + u>i + u 2 dk (2u> k + Lu 2 )(2u> k — u> 3 ) 

d_ p(k) d (3(k) 

dk (2oj k + uJi)(2uj k — uo 2 — uj 3 ) dk 2u; k — uo 2 

_d_ (3{k) d (3(k) 

dk (2uj k — aj 3 )(2c<jfc + lu\ + uj 2 ) dk 2uj k + uj\ 

f _a m d m 

dk (2u> k + u>3)(2uj k — loi — uj 2 ) dk 2uj k — uj\ 

d_ (3{k) d j3(k) 

dk (2ui k — u! 3 )(2ujk + u>i + lo 2 ) dk 2u k + u>i 



+ 



(B19) 
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Since the above expression is an even function of u>k , we obtain 



h + h + l7 = 0- (B20) 



2. 



3. Let 



Then 



u Uk) = i d m 

2 dk (2ujk — Wi — w 2 )(2tjfe — lu± — LU2 — W3) <9fc 2wfc — wi 

_ /3(fc)_a_ 1 _a_ ^(fc) 

2 dk 2u>k — u>\ — u>2 — ^3 <9fc (2cijfc — uj-\)(2u>k — u>\ — W2) 
/3 2 (fc) r _g 1 g 1 



2wfc — u>i dk lu — lui — lu 2 + uik dk lu — lui — lu 2 — ^3 + 



2uik — u)\ K dk u — u)\ — u>2 + ujk dk ui — ui\ — lo 2 — 0J3 + u>k ^ k 
d 0(k) d 1 



+ P(k) 



+ (3(k) 



_ dm k /3 2 (k) d 1 

dk 2(2ujk — u>i) 2 dk (2uik — u>± — u>2){2ujk — u)\ — L02 — ^3) 
du k f3 2 {k) 1 d 1 

dk 2(2uJk — t^i) 2 {2uok — Wi — ^2) dk 2u>k — uji — uj 2 — u>3 

_ du k (3 2 (k) 1 d 1 

dk 2(2uJk — oj\) (2ujk — — UJ2) 2 dk (2iVk — uj\ — 002 — W3) 

_ duj k (3 2 {k) 1 d 1 

dk (2bJk - lui) {2uj k — wi - w 2 )(2o) fe - ui - u>2 — W3) dk 2uj k - u\ - lu 2 

:= Ti + T 2 + T 3 + T 4 + T 5 + T 6 + T r + T s + T 9 , 



(B21) 



T{uj k ) = I 2 {u) k ) + J 3 (w fc ) + / 4 (w fc ). (B22) 

_(3(k)_d_ 1 d p{k) 

2 dk (2uk — oji — uj 2 )(2ujk — u)\ — 002 — ^3) dk 2u>k — Wi 

(3(k) d 1 d_ f3(k) 

2 dk 2u>k — <*>i — u>2 — ^3 dk (2u>k — ^>i){2u>k — u>i — L02) 
2(3 2 {k) d_ 1 d_ 1 



dk (u> — u>i + LOk)(^ — u>i — u) 2 + Wfe) dk u> — u>i — LU2 — LU3 + LUk ^ U>k 

d 13(h) d 1 J 

dk uj - uii + ujk dk (w -UJ1-0J2 + u>k)(u> - uj\ - lo 2 - ^3 + Vk) U ~^ h (B23) 

?2 



where 
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T 3 + T 4 + T 5 = +0{k) 



+0(k) 



( d 0(h) \ ( 1 d 



Then 



\ 8k U — U>1 + U>k J \LU — LUl — UJ2 + dk U> — LUl — LU2 — LU3 + LL>k , 

d 0(k) \ ( 8 1 



\U)=UJ k 



= + 
+ 



dk lo — u>i + LJk J \dk (uj — u>i — u>2 + ) (w — u>\ — u>2 — L03 + Wfc) 
2 {k) d 1 8 1 

2u>k — oji dk u — u)\ — u>2 — 0J3 + uok dk uo — uj\ — L02 + ujk u} ~ ulk 

d0(k) 0{k) d_ 1 

dk 2(2ujk — uj\){2uLik — Wi — 0J2) dk 2uk — — ^2 — ^3 

d0(k) 0(k) d 1 

dk 2(2cjfe — u>i) dk (2ujk — u>i — U2)(2ujk — u>i — L02 — W3) 

duj k 2 {k) 1 d 1 

dk 2(2cjfc — ujx) 2 (2uuk — coi — W2) dk 2u>k — u>\ — 0J2 — ^3 

duj k 2 {k) d 1 

dk 2(2u>k — ui) 2 dk (2u>k — u>\ — u>2)(2ujk — u>\ — L02 — W3) 
2 {k) ( d 1 \( 8 1 



4(2wfe — uji) \dk 2u>k — oji — ^2 ) \dk 2uk — loi — L02 — ^3 

2 {k) 1 d 2 1 

2(2ujk — 0J1) {2uok — u\ — U02 — ^3) dk 2 2tdk — lu\ — LU2 

doj k 2 {k) 1 d 1 

dk (2u>k — 0J\) (2u!k — wi — L02){2u>k — wi — ll>2 — UJ3) dk 2u>k — u>i — u>2 



T' = T 3 + T 4 + T 5 +T 6 + T 7 + T 8 + T 9 

80(k) 0{k) d 



= + 
+ 
+ 
+ 
+ 
+ 

= + 
+ 
+ 



dk 2(2cjfe — wi)(2wfc — uji — U02) dk 2tOk — wi — ^2 — ^3 
d0(k) 0{k) d 1 

dk 2{2tOk — uoi) dk {2uok — — uo2){2ujk — uji — u)2 — 003) 
2 (k) ( 8 I \ ( 1 d 1 



2 \dk 2ujk — ui\) \ 2u)k — 0J1 — u>2 dk 2u>k — lo\ — u>2 — 0J3 

2 {k) ( d_ 1 \ ( d_ 1 

2 \dk 2tOk — wi J \dk (2u>k — wi — ^2){2ojk — <jJi — 1^2 — ^3) 
2 {k) ( d 1 \(d 1 



2(2cjfe — cji) \dk 2ujk — u>\ — u>2 ) \8k 2uJk — uji — 0J2 — ^3 

2 {k) 1 d 2 1 

2(2cjfe — uj\) (2u>k — u>i — u>2 — UJ3) dk 2 2ujk — u\ — u>2 
0{k) ( d 0{k) \ ( 1 d 1 



2 \dk 2u!k — ui) \2u>fc — oji — u>2 dk 2uik — ^1 — ^2 — ^3 

m (9 0{k) \ ?d_ 1 

2 \dk 2u>k — w\) \dk {2uok — ^1 — ^2){2uok — uoi — 002 — ^3) 

2 {k) d 1 d 1 

2(2cjfe — uoi) dk 2cdk — u\ — u 2 — 003 dk 2tOk — uo\ — u 2 



lLU = UJ k 



Therefore 
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T(w fc ) = T 1 +T 2 + T / 

(3(k) 1 d 2 p(k) 

2 (2u> k — u>i — u) 2 ){2u> k — loi — ll> 2 — W3) <9fc 2 2u^ — wi 

g(fe) 1 j¥ P(k) 

2 2uj k - ui - w 2 - dfc 2 (2^fc - wi)(2wfc - u>i - u> 2 ) 

(P(k) d 2 1 

2(2o>fe — wi)(2wfe — wi — c^2 — W3) 9fc 2 2wfe — wi — 

/3(fc) d 2 I3{k) 

(2uj k - uji - uj 2 )(2u k -W1-W2- W3) <9/c 2 2w fc - Wi 
/3(fc) / d f3(k) \ ( d 1 



+ 



2oJk — u>\ — u>2 — d>3 \dk 2ujk — u>i ) \dk 2uj k — uj\ — u 2 

f3(k) d 1 d (3{k) 

2uj k — Wi — uj 2 — dk 2ui k — u>i — oj 2 dk 2uj k — Wi ' 



Combining the above results, we obtain 



l'(u)i,W2,uz) = h(u)k) + T(u)k) - h(-Uk) - T(-u k ) 

/3 4 (fc)(4co fe - Wl -^ 3 ) 



(2uj k - u>i)(2u>k + uj 2 )(2uJk - ^3)(2w fe - wi - u>2 - W3) 
P\k){Auj k +lo 1 + u 3 ) 

(2w fe + wi)(2w/s - W2)(2o;fe + uj 3 )(2uj k + lji + u 2 + L03) 
(3{k) d 1 d pjk) 

2u> k — loi — lo 2 — LU3 dk 2uj k — ui\ — uj 2 dk 2ui k — u\ 
(5{k) d 1 d (3{k) 

2ui k + uu\ + uo 2 + dk 2uj k + 0^1+^2 dk 2ui k + u)\ ' 

and 

S=l ]T 2"(wi,W2,ws) (B24) 

V(u>\,U>2,U>3) 

We now proceed to compute the general four-wave-mixing third order susceptibility X^ 3 H u; i> a; 2j w 3) which is defined 
by: 

x (3)( Wl , W2 , a , 3 ) = H^ol^;5. (B25) 
For infinite chains, the summation is replaced by the integral: 

D 4„, Ah 0»4, 



(3), , 2e 4 n dk 2e 4 n 



We first recall: 



1 2e 4 n /"^ 

6 £ J Z'("u"2,u 3 )dk (B26) 

X(W1,W2,W 3 ). 



= (B27) 



e(k) = y/[2t cos (fca)] 2 + [A sin (fca)] 2 , (B28) 
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«*>--?& 

Furthermore, we define the following quantities: 

6 ~ 2h' 



(3) _ 8 e 4 n a 3 



Let us introduce change of variable: 



Then for < k < we have: 



and hence 



Zl = ~2~X' fori = 1 ' 2 ' 3 ' 



{ 



_ q ^(l-<5 2 x 2 )(x 2 -l) 
dk 5 x 



1/5 xdx 



y/(l-S 2 X 2 )(x 2 -l) 

(2x - Zl - z 3 ) 



x 8 (x - zi)(x + z 2 )(x - z 3 )(x - zi - z 2 - z 3 ) 
(2x + zx+ z 3 ) 



x 8 (x + z\){x - z 2 ){x + z 3 )(x + zi + z 2 + z 3 ) 

+ 



2S\ 2 1 9 19 1 



a J x 2 (x — zi — z 2 — z 3 ) dk x — z\ — z 2 dk x 2 (x — z\) 

25\ 2 1 9 19 1 

a J x 2 (x + z\ + z 2 + z 3 ) dk x + zi + z 2 dk x 2 (x + z\) 

(3) 15 f 1/s xdx 



X ° 1024 J i ^(l-5 2 x 2 ){x 2 -1) 

r (2x - zi - z 3 ) 

x s (x - zi)(x + z 2 )(x - z 3 )(x - zi - z 2 - z 3 ) 

(2x + zi + z 3 ) 
x 8 (x + zi)(x — z 2 )(x + z 3 )(x + zi + z 2 + z 3 ) 

+ 



2S\ 2 1 / 9 1 \ / 9 1 



a J x — zi — z 2 \dk x 2 (x — z\) ) \dk x 2 (x — zi — z 2 — z 3 ) 



+ — 



2S\ 2 1 / 9 1 \ / d 1 



a J x + zi + z 2 \dk x 2 (x + z\) ) \dk x 2 (x + zi + z 2 + z 3 ) 
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By (B34), we have 



(3) 15 f 1/s xdx 



1024 J 1 ^(l-6 2 x 2 )(x 2 -1) 

(2x - zi - z 3 ) 

x 8 (x - zi)(x + z 2 ){x - z 3 )(x - z\ - z 2 - z 3 ) 

(2x + gi + z 3 ) 

x a (x + z\)(x - z 2 ){x + z 3 )(x + z\ + z 2 + z 3 ) 
A{l-5 2 x 2 ){x 2 -1) / d 1 \ ( d 1 



x 2 (x — z\ — z 2 ) \dx x 2 (x — zi) J \ dx x 2 (x — z\ — z 2 — z 3 ) 
4(1 - S 2 x 2 )(x 2 - 1) ( d I \ ( d 1 



(3) 

= Xo 



x 2 (x + z\ + z 2 ) \dx x 2 (x + z\) J \dx x 2 (x + z\ + z 2 + z 3 ) 

L 



15 r 1/s xdx 

1024 J 1 ^(l-S 2 x 2 )(x T ^T) 

(2x - zi - z 3 ) 

x 8 (x - zi){x + z 2 ){x - z 3 )(x - zi - z 2 - z 3 ) 

(2s + zi + z 3 ) 

x 8 (x + zi)(x - z 2 ){x + z 3 )(x + z\ + z 2 + z 3 ) 

4(1 - 5 2 x 2 ){x 2 - 1) (3a; - 2zi)(3x - 2{ Zl + z 2 + z 3 )) 



x 8 (x — Z\ — z 2 ) (x — Z\) 2 (x — Z\ — z 2 — z 3 ) 2 
4(1 - S 2 x 2 )(x 2 - 1) (3x + 2zx)(3x + 2{z x + z 2 + z 3 )) \ 
x 8 (x + z 1 + z 2 ) (x + z 1 ) 2 (x + Z\ + z 2 + z 3 ) 2 J 
:= h(uj±,u>2,u> 3 ) + J 2 (wi, u>2, u 3 ) + I 3 {u>i,uj 2 ,uj 3 ) + J 4 (wi, w 2 , w 3 ), 



With the aid of Mathematica, we find that 



J= (Ii(uji,uj 2 ,uj 3 ) + I 2 (uji,uj 2 ,uj 3 )) 

15x^ 3) f 1/S dx 



1; 



128(«i + z 2 ){z 2 + z 3 ){z 3 + *i) J 1 ^(1- 5 2 x 2 )(x 2 -1) 
_j_/y- z[_ (z 1 +z 2 + z 3 f 

T 8 \ Z-^i T 2 _ 



X 2 - Z 2 X 2 - (Zl + Z 2 + Z 3 ) 2 



and 
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15x 3) f 1/S y/{l-5 2 x 2 ){x 2 - l)dx 



128 J 1 x s 

E (zi + z 2 ) 5 (z 1 + z 2 - 2z 3 ) 
... =l zl44(x* + 

-{Zl +Z 2 + Z 3 ) 3 ((Z1 + Z 2 + Z 3 f{ Zl Z 2 + Z 2 Z 3 + g3£l) - z 1 z 2 z 3 (8(z 1 +z 2 + z 3 ) 3 + 7ziz 2 z 3 )) 

z\z\zl{zx + z 2 ){z 2 + z 3 ){z 3 + Zl )(x 2 - ( Zl + z 2 + z 3 ) 2 ) 
(zi + z 2 + z 3 ) 3 ((zi + z 2 + z 3 f + z x z 2 z 3 )(x 2 + (zi + z 2 + z 3 ) 2 ) 



+ 



z 1 z 2 z 3 (z 1 + z 2 )(z 2 + z 3 )(z 3 + zx)(x 2 - (z x + z 2 + z 3 ) 2 ) 2 

y, -Z 2 ( Zl - 2(Z 2 + Z 3 ))(X 2 + Z\) 

z 2 z 3 (z 2 + z 3 )(x 2 - z 2 ) 2 



E 



-z\{z\{z 2 + z 3 ) - 8z 2 z 3 (z 2 + z 3 ) - 2zi(z 2 + z 3 ) 2 + lz x z 2 z 3 ) 



z 2 z 2 {z 2 + z 3 ){x 2 - z 2 ) 
Using integration by parts, we obtain 

" 1/S ^(1-S 2 x 2 )(x 2 -l)dx x 2 + z 2 



J l x s x 2 -z 2 V (l-8 2 x 2 ){x 2 - 1) 

f _ 



^(l-^x 2 )^ 2 -!)^ 1 ( 7 | 1 - ,5 2 x 4 



(1 - 5 2 x 2 )(a; 2 - 1) 



Therefore, 



15x[, 3) /- 17,5 ^(l-^x 2 )^ 2 -!)^ 



if = 



128 Jj x s 

" (z 1 + z 2 ) 5 (z 1 + z 2 ~2z 3 ) 
~i z\zlzl(x 2 - (Z! + z 2 ) 2 ) 

-{ Zl + z 2 + z 3 f 



E 



+ 



+E 



z\z\z\(x 2 - (Z! + z 2 + z 3 ) 2 ) 

(£i + z 2 + z 3 ) 3 ((zi + z 2 + z 3 ) 3 + Zxz 2 z 3 ) 1 - <5 2 x 4 

z 1 z 2 z 3 (z 1 + z 2 )(z 2 + z 3 )(z 3 + Zl) (x 2 - {zi + z 2 + z 3 ) 2 )(l - S 2 x 2 )(x 2 - 1) 

■z 2 ( Zl - 2(z 2 + z 3 )) 1 - 5 2 x 4 



^ z 2 z 3 (z 2 + z 3 )(x 2 - z\) (1 - <5 2 x 2 )(x 2 - 1) 

E -z\{z\ - 2zi(z 2 + z 3 ) + 6z 2 z 3 ) 
z 2 2 z 2 {x 2 ~z 2 ) 



(B38) 



(x 2 — z 2 ) 2 

- '* y/(l-6*x*)(x*-l)dx 1 ( ^x^l + S 2 -^)' {m)) 



(B40) 
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Finally, 

X (3) = J + K 



15 X ^ dx 



L 



128 J 1 x»^(l-S 2 x 2 ){x 2 -1) 
Z 3 z 3 -(Z 3 + z 1 z 2 z 3 )5 2 ^ z 2 {-z 1 + 2{z 2 + z 3 )){l-5 2 x i ) 



+ 



z\z 2 z 3 o x 2 -Z 2 j-^(j{x 2 - z 2 ) p(z f^ 2 Z3) 2z 2 z 3 (z 2 + z 3 ){x 2 - z 2 ) j ( B41 ) 

15 X o } f 1/S V(l~5 2 x 2 ){x 2 -l)dx 



i; 



128 J 1 x 8 

y- / (z 1 +z 2 ) 5 (z 1 +z 2 -2z 3 ) -z 2 (z 2 - 2z x {z 2 + z 3 ) + 6z 2 z 3 ) \ Z 6 
^ \ 2z\z 2 z 2 {x 2 - ( Zl + z 2 ) 2 ) 2z 2 z 2 (x 2 - z\) J 



.H:,. --'-■>■ ' - - ! ' ' Z f Z 2 Z U x2 Z2 ) 

where 

a := ( Zl + z 2 )(z 2 + z 3 )(z 3 + Zl ), (B42) 

and 

Z := zi + z 2 + z 3 . (B43) 

We now consider the case 5 — > 0. 

Using the results in Appendix A, we obtain 



.(3) _ LO M) 



12© I ZiZ 2 Z 3 <7 z^z 2 z 3 z — j a 

2/ 



+ e zi( r i+ , 2(za+ f 3)) ^(4,^) 

P(™) 2 ^2 + . 3 ) 

, 15Xo 3) / V- (z 1+ z 2 ) 5 (z 1+ z 2 -2z 3 ) 
+ ^ 2. 2^ M(4,z 1 + z 2 ) 

-2^2 



-z\{z\ - 2z 1 ( z 2 + z 3 ) + 6z 2 z 3 ) 



(B44) 
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FIG. 1. The Feynman diagram for x^ 3 ' (— ^, — K; ui\, ki, LU2, k 2 ,o;3, k 3 ), where Q = lo\ + L02 + <^3, K = ki + k 2 + k 3 . 
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FIG. 2. (Color) The magnitude of four-wave-mixing(FMW) x' 3 ' ^li ^2, W3) under SSH model is in unit of 10 8 esu. Z% 
is defined by Eq.(3.5). 
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FIG. 3. The hyperpolarizabilities under the SSH model for frequency region < z < 2.0, where z = hu>/(2A): DCKerr (top 
left), DCSHG (top right), IDIR (bottom left) and EFIOR (bottom right). 
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FIG. 4. The hyperpolarizabilities under the SSH model for non-resonant region: 
IDIR (long dashed line), DCSHG (dashed line) and THG (dot-dashed line). 
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FIG. 5. The hyperpolarizabilities under the TLM model for non-resonant region: Xtlm(~ w i w i w ) (real line) versus 
Xtlm( uj > lu > uj > —3a;) (long dashed line), where Z = huj/2A. 
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